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The apparent anomaly in the ratio of muon to electron atmospheric neutrinos first observed 
by Kamiokande and 1MB has been confirmed by Super-Kamiokande and Soudan-2. The 
experimental analysis, including the asymmetry in the zenithal distributions of the /i — type 
events in Super-Kamiokande gives a strong support to the neutrino oscillation hypothesis to 
solve the anomaly. 

In this work we are interested by the role of nuclear physics in the neutrino-oxygen reactions 
used to detect the atmospheric neutrinos. We point out that multi-nucleon excitations of 
np-nh type and that nuclear correlations could modify an experimental analysis a la Super- 
Kamiokande because they lead to a substantial enhancement of the number of 1 Cerenkov 
ring retained events. 
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1 Introduction 

The apparent anomaly in the numbers of f^ and z^e prgicj^^d by interactions of cosmic rays in the 



atmosphere has been confirmed by several experiments liloU'Q. While the expected ratio of u^ to Ve 
should be of the order of two in favour of the muonic neutrinos, as indicated bv simple arguments on 
the decays of pions and kaons and confirmed by several theoretical calculationsSfl, the measured ratio 
is significantly lower. Indeed the double ratio: 
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between the experimental and the simulated u^ to Ve ratios (R^/e) is found to be of the order of ~ 0.6. 
Furthermore the Super-Kamiokande collaboration has observed an asymmetry in the zenithal distri- 
butions of the z/^ events (both sub- and multi- GeV) 0. It appears that the number of upward-going 
Ufj^ (i.e. coming from the antipodes) is smaller than expected while the number of downward-going v^ 
corresponds to the expectations. Together with the symmetric behaviour of the zenithal distributions 
of the z/g events and the results of the Chooz experiment which excludes the possible f^ — > Ue solution 
to the Super-Kamiokande results, this favours a solution of the atmospheric neutrino problem in the 
f^ — > ux oscillations channels, with vx = z^r or ux = I's (some kind of sterile neutrino). 



2 Role of nuclear physics 

In this work we are interested in the experiments using large underground water Cerenkov detectors: 
Kamiokande and Super-Kamiokande, 1MB. The analysis is based on events where only one Cerenkov 
ring is detected. These "one Cerenkov ring" (1 C.R.) events are usually assumed to be produced 
by quasi-elastic charged current interactions in which a charged lepton is emitted above Cerenkov 
threshold and leads to the Cerenkov ring. The nucleon which is ejected from the nucleus is in general 
below threshold and therefore does not produce another ring (in water the threshold kinetic energy for 
a nucleon is about 0.5 GeV and the assumption is rather correct). However we know that the nuclear 
dynamics is far more complex than this simple picture and there are others sources of 1 C.R. eventsQ. 
Indeed the region of energy transfer in processes involving atmospheric neutrinos of ~ 1 GeV extends 
from the quasi-elastic peak to the Delta resonance region. The quasi-elastic region provides the most 
important part of the 1 G.R. events with the previous assumption on the emitted nucleon. At higher 
energy we know from electron scattering that 2 particles - 2 holes excitations are necessary to fill the 
intermediate region (the so-called dip region) between the quasi-elastic and the Delta resonance peaks. 
In the context of atmospheric neutrino such excitations lead to the emission of two nucleons and one 
charged lepton. Assuming that the nucleons are under threshold one still gets an 1 C.R. event. The 
same arguments are valid for n particles - n holes excitations. 

At still higher energy we evaluate the nuclear responses in the Delta resonance region. If the Delta 
decays into a pion and a nucleon and if the pion escapes from the nuclear medium above threshold 
('^ 70 MeV in water for a pion), then one gets two Cerenkov rings (one for the charged lepton and 
the second for the emitted pion). Such events are rejected by the experimental cuts. However the 
pion in the nucleus is a quasi-particle with a broad width and can decay for instance into a particle- 
hole excitation. Therefore the decay of a Delta in the nuclear medium can lead to a nucleon and a 
nucleon-hole state. Such a process, where two nucleons are ejected from the nucleus, just produces 
one Cerenkov ring. 

The previous arguments lead to the conclusion that a full calculation of the neutrino-oxygen cross 
sections beyond the quasi-elastic assumption with identification of the possible final states is necessary. 
This procedure gives the total 1 C.R. events yields in the atmospheric neutrinos experiments which 
have to be compared to the quasi-elastic 1 C.R. events yields usually retained in the simulations. Note 
that our investigation relies on very simplifying assumptions concerning the particles in the final state. 
In particular we will neglect throughout this paper the role of pion re-interactions and absorption, the 



emission of nucleons above Cerenkov threshold and the possible particle misidentification problems. 
This last assumption is due to a lack of information on the experimental detection efhcencies. 

3 Nuclear responses and differential cross sections 

The starting point of this calculation is the inclusive charged-current cross section for the reaction: 
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where Gi;' is the weak coupling constant, ^c the Cabbibo angle, k and A;' the initial and final lepton 
momenta, q^ = k^ — k' = (oj, q) the four momentum transferred to the nucleus, 9 the scattering 
angle. Ma (M) the Delta (nucleon) mass. The plus (minus) sign in eq. (|2|) stands for the neutrino 
(anti-neutrino) case. In a provisional approximation we have neglected in eq. (|2|) the lepton masses 
and we have kept the leading terms in the development of the hadronic current in p/M, where p 
denotes the initial nucleon momentum. The electric, magnetic and axial form factors are taken in 
the standard dipole parameterization. We have introduced the inclusive isospin (Rr), spin-isospin 
longitudinal {RaT{L)) ^"^^ spin-isospin transverse {R^^^rp^) nuclear responses functions corresponding 
to the isospin operator (tq) and to the different projections of the spin-isospin operators (cr r^) on 
the direction of the momentum transfer (q). Assuming that 0^{j) is one of the previous operator, 
responsible for the coupling of the incident neutrino to a Particle-hole state, the corresponding nuclear 
response is given by: 

Ra"" = E (-1 E O^ij) e^^-^^ |0) {n\ J: O^'ik) e^9-^'=|0)*5(a; -E^ + Eo) (3) 

n j=l k=l 

The challenging task is the evaluation of the various nuclear responses. The main problem arising is 
the knowledge of the nuclear excited states for transfer energies of several hundred MeV. We avoid 
this difficulty with the model developped by Delorme and GuichonB which is based on a semi-classical 
approximation. This implies the use of local Fermi momentum kpir) which is calculated by the means 
of an experimental nuclear density: kp{r) = (3/27r^ p(r))^'^ and this leads to the "bare" nuclear 
responses I^{uj,\q\) (in the following "bare" will mean that the nuclear correlations are switched off). 
Indeed until now we have considered that the particle-hole excitations in the medium do not interact 
with each other. This is of course not true. To properly take this effect into account we adopt a 
phenomenological parameterization of the interaction between excitations in each reaction channel: a 
pure contact force with a Landau-Migdal parameter /' is taken in the isospin channel and a (vr -|- p)- 
exchange plus a contact force with a Landau-Migdal parameter g' in the spin-isospin channel. Once 
the form of this potential is known, we exactly solve the RPA equations in the ring approximation. 
This gives the full RPA nuclear responses i?(a;, \q\). 

The modified Delta width in the nuclear medium is split into the contributions of different decay 
channels 13: the "quasi-elastic" channel, A — > ttN, modified by the Pauli blocking of the nucleon 
and the distorsion of the pion, the two-body {2p-2h) and three-body {3p-3h) absorption channels. 
The {2p-2h) excitations which are not already included in the modified Delta width are evaluated by 
extrapolating the calculations of two-body pion absorption at threshold 113. These parameterizations 
lead to a good description of the pion-nucleus reactions. 



The inclusive nuclear response functions are split into their different contributions. By construction, 
the "bare" response is the sum of the quasi-elastic, the n particles - n holes (n = 2,3) and the resonant 
Delta (A — > n N) responses. In the RPA case we will recognize the same partial contributions 
supplemented by the coherent pion emission channel which arises when an intermediate pion exchanged 
between excited states is placed on its mass-shell. The full response can then be written as: 



R{UJ, q) = Rqel{uJ, q) + Rnp-nhW, q) + Rresi^, q) + Rcoh{^, 



(4) 



We recall that within our model Rqc,i and Rnp-nh will induce 1 C.R. events and that Rres and Rcoh 
will lead to the emission of 2 Cerenkov rings. 

The results for the differential cross section da/dk' , which is obtained from the doubly differential 
cross section (eq. (|2|)) by a numerical integration over the solid angle, are shown on figure (|l|). The 
neutrino energy is fixed at the typical value E^ = 1 GeV. 
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Figure 1: Differential charged-current Ue —^^ O interactions cross-section versus the energy transfer. The thick curve 
represents the inclusive cross-section. The following exclusive contributions to the inclusive cross-section are displayed: 
quasi-elastic (full thin curve), (np-nh) (dotted curve), resonant Delta (A — > nN) (short dashed curve). Also shown are 
the "bare" (without RPA correlations) quasi-elastic (long dashed curve) and the coherent pion emission cross-sections 

(dot-dashed curve). 

The inclusive cross section is given by the thick curve. It gets its main contribution from the quasi- 
elastic channel (thin full line) which peaks at relatively low energy transfer. For a sake of comparison 
we have shown the contribution of the quasi-elastic channel without RPA (thin long dashed line) . We 
observe that the cross section is reduced and hardened in the RPA case. The shift in strength is due 
to the repulsive RPA-interaction in the spin-isospin transverse response. Indeed the Landau-Migdal 
interaction is repulsive for all values of the transfer and the p — exchange piece is not attractive enough 
in the domain of energy considered here to counteract this feature. For instance we know that this 
conclusion no longer holds in the longitudinal channel where the vr — exchange is attractive enough to 
overcome the Landau-Migdal interaction and to create a collective mode (the pionic brancht3). But, 
as we already see from eq. (§), the spin-isospin longitudinal channel is suppressed with regard to the 
transverse one and it plays a minor role in the neutrino-nucleus reactions. 



The effect of the RPA correlations are less strong in the others partial channels and are not shown 
here. The resonant Delta channel (short dashed curve) arises at high energy transfer {oo ~ 450 MeV). 
The most intereting feature of the cross section is the importance of the (np-nh) channel (dotted curve). 
We see that it gives a rather large contribution to the inclusive cross section and has an extended 
spectrum from the quasi-elastic to the Delta peaks. In particular these processes are important in the 
dip region, as mentionned previously. 

The conclusion is that the inclusive neutrino-oxygen cross section is strongly modified with respect to 
the "bare" quasi-elastic one, which is often the sole channel entering into the simulations. In particular 
we have seen the occurence of large contributions from the two- and three-body channels. 

4 Neutrino- Oxygen events yields 

The interesting quantities are the neutrino-oxygen events yields, which we compute at fixed charged 
lepton momentum: 

Yi.^ + f?^){k') = Q dE {^^.^E)^-^^^^^^{E,k') + ^^AE)^^^^Q^{E,k')^ (5) 

where E is the neutrino energy, ^i, the incoming neutrinos flux and da/dk' the neutrino-oxygen cross 
section computed in the previous section. We use the fluxes of Bartol in our calculations. The 
main feature of these fluxes is their sharp decrease with increasing neutrino energy. Within the few 
rough assumptions described previously, we classify our events with the number of Cerenkov rings 
they produce. The results for the 1 C.R. events yields, which are relevant in the atmospheric neutrino 
experiments, are shown on fig. (^) for incident i/^ and i^^, the full curves corresponding to the total 1 
C.R. events yields and the dashed curves to the sole quasi-elastic 1 C.R. events yields. We give the 
results of the calculations without (thin curves) and with (thick curves) RPA.f] 

First we observe that the RPA tends to reduce the events yields. Indeed it shifts the cross section 
towards high energies which are disfavored by the neutrino flux. The maximal reduction factor is of 
the order of 10 %. As expected, the {np-nh) excitations on the contrary increase the absolute events 
yields. At the maximum value of the yields, the enhancement of the total yield with respect to the 
"bare" quasi-elastic one is around 20 %. This result reflects the main features of the cross sections. 
Furthermore one must be aware that this result is just a lower limit of the true enhancement. Indeed 
if a pion produced through Delta decay or coherent emission is absorbed in the medium, it will not 
produce a Cerenkov ring and therefore the whole corresponding charged-current process belongs to the 
1 C.R. events class. This problem remains to be investigated. Nevertheless we can conclude that the 
neglected reactions of the atmospheric neutrinos on pairs of correlated particles in an oxygen nucleus, 
which lead to 1 C.R. events, are a possible explanation for the under-estimation of the experimental 
^e + ^e data by the Monte-Carlo simulations in Super-Kamiokande |. Furthermore we insist on the 
fact that these modifications to the theoretical events rates could have rather large consequences on 
the results of the experimental collaborations which are now based on the absolute rates of events and 
not only on the ratios of these events rates. 



5 Conclusion 

In this work we studied the neutrino-oxygen reactions in the domain of energy relevant for the at- 
mospheric neutrinos experiments using large underground Cerenkov detectors. First we achieved the 
calculations of the nuclear responses of oxygen to a neutrino excitation. The semi-classical model used 



Let us note that in order to give a full calculation of the events rates, one has to include the neutral currents. In 

this case, because the scattered neutrino is not seen experimentally, the resonant and coherent channels belong to the 

required 1 C.R. class (assuming that the pion is above threshold). The calculations have been performed but the analysis 

needs more informations on detection efhcencies and on absorption of pious. This problem has to be investigated further. 

"^We thank Pr. Y. Declais who caught our attention on this point. 
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Figure 2: One Cerenkov r^^ + O events yields versus the muon momentum. The full curves correspond to the total 1 

C.R. events yields with (thick curve) and without (thin curve) RPA. The dashed curves correspond to the quasi-elastic 

1 C.R. events yields with (thick curve) and without (thin curve) RPA. 

here has been developped to study the nuclear intermediate energy spectrum, from the quasi-elastic 
to the Delta peaks. It includes a parameterization of the multi-particles excitations {np-nh, n = 2,3) 
and the possible interactions between excited states by the means of the ring approximation to the 
RPA. 

Within a few assumptions we perform an analysis of our results with respect to the number of Cerenkov 
rings produced in the various nuclear processes considered. We conclude first that the RPA correlations 
tend to reduce the quasi-elastic and therefore the total 1 Cerenkov ring events yields, the maximal 
reduction factor being of the order --^ 10 %. These events are those retained by the experimental cuts 
performed in the Kamiokande, 1MB and Super-Kamiokande collaborations, whose results are our main 
point of interest. The second result is that the {np-nh) excitations increase the total 1 Cerenkov ring 
events yields with respect to the "bare" quasi-elastic one usually taken into account. The strong 
enhancement factor (~ 20 %) could lead to modifications in analyses of experimental results based on 
absolute number of events. 
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Appendix 

Effects of the RPA on the differential (figure |3|) and total (figure §) i'^{i'^) +^^ O cross sections. The 
total cross sections are obtained by numerical integration of the differential ones da/dk' on the lepton 
momentum. The thick curves display the RPA cross sections while the thin dashed curves represent 
the bare cross sections. The Vfj, +^^ O are the most sensitive reaction channels to the RPA effects. 
This is true for the whole range of neutrino energy considered here and results in a quenching of the 
cross sections which is less than 10 %. 

The last figure (figure ^) shows the various contributions to the inclusive total cross sections. The 
RPA responses ar used to evaluate the cross sections. The most important cross sections are the 
quasi-elastic and the resonant Delta ones (respectively full thin and dashed curves). The rather large 
contribution of the np — nh reaction channels (dotted curve) is obviously seen on the figure. Finally 
the coherent pion emission channel (dot-dashed curve) tends to increase with the neutrino energy. 
However one has to notice that the contribution of this channel remains rather weak (typically one 
order of magnitude below the previous contributions). 



v - O differential cross sections (E = 1 GeV) 



Total RPA 
Total bate 




600 700 

CO (MeV) 



Figure 3: !^^(!^^)+^®0 differential cross sections versus the energy transfer. The full thick (dashed thin) curves correspond 

to the RPA (bare) calculations. 
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Figure 4; !^^(i^^) +^^ O total cross sections versus the neutrino energy. The full thick (dashed thin) curves correspond to 

the RPA (bare) calculations. 
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Figure 5: u^^ +^^ O total RPA cross sections versus the neutrino energy. The full thick curve correspond to the inclusive 

reaction channel which is the sum of the following partial channels: quasi-elastic (full thin curve), Delta resonant (dashed 

curve), np — nh (dotted curve), coherent pion emission (dot-dashed curve). 



